Band-gap solitons in nonlinear optically-induced lattices 
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We introduce novel optical solitons that consist of a periodic and a spatially localized components 
coupled nonlinearly via cross-phase modulation. The spatially localized optical field can be treated 
as a gap soliton supported by the optically-induced nonlinear grating. We find different types of 
these band-gap composite solitons and demonstrate their dynamical stability. 



Recent theoretical and experimental results demon- 
strated nonlinear localization of light in optically-induced 
refractive index gratings 0, Q] ■ Such localized states can 
be treated as "discrete" and "gap" solitons observed in 
fabricated periodic photonic structures , but supported 
by gratings induced by a complementary optical field. 
Optically induced lattices open up an exciting possibility 
for creating dynamically reconfigurable photonic struc- 
tures in bulk nonlinear media. The physics of coherent 
light propagating in such structures can be linked to the 
phenomena exhibited by coherent matter waves (Bose- 
Einstein condensates) in optical lattices 0. 

Among the most challenging problems in the physics 
of induced gratings is the creation of stable, uniform pe- 
riodic optical patterns which can effectively modulate the 
refractive index of a nonlinear medium. Periodic modu- 
lation of the refractive index can be induced, for instance, 
by an interference pattern illuminating a photorefractive 
crystal with a strong electro-optic anisotropy Inter- 
fering plane waves modulate the space-charge field in the 
crystal, which relates to the refractive index via electro- 
optic coefficients. The latter are substantially different 
for the two orthogonal polarizations. As a result, the ma- 
terial nonlinearity experienced by waves polarized in the 
direction of the c - axis of the crystal is up to two orders of 
magnitude larger than that experienced by the orthogo- 
nally polarized ones. When the lattice-forming waves are 
polarized orthogonally to the c - axis, the nonlinear self- 
action as well as any cross-action from the co-propagating 
probe beam can be neglected. The periodic interference 
pattern propagates in the diffraction-free linear regime, 
thus creating a stationary refractive-index grating [2| • 

In this Letter, we develop the concept of optically- 
induced gratings beyond the limit of weak material non- 
linearity and propose the idea of robust nonlinearity- 
assisted optical lattices, created by nonlinear periodic 
waves. Strong incoherent interaction of such a grating 
with a probe beam, through the nonlinear cross-phase- 
modulation (XPM) effect, facilitates the formation of a 
novel type of a composite optical soliton, where one of 
the components creates a periodic photonic structure, 
while the other component experiences Bragg reflection 
from this structure and can form gap solitons localized 



in the transmission gaps of the linear spectrum. The 
observation of nonlinear light localization in this type 
of optically-induced gratings can be achieved in photore- 
fractive medium with two incoherently interacting beams 
of the same polarization. We study such a configuration 
in a saturable medium and demonstrate the existence and 
stable dynamics of these novel band-gap lattice solitons. 

The propagation of two incoherently interacting beams 
in a photorefractive crystal can be approximately de- 
scribed by the coupled nonlinear Schrodinger (NLS) 
equations for the slowly varying envelopes £"„ {n = 1,2), 



dz 



aN{I)Er, = 0, 



(1) 



where N{I) — I /{I -I- si) describes saturable nonlinear- 
ity, I = J2 l^nP is the total light intensity, s is the sat- 
uration parameter, and a = ±1 stands for the focusing 
or defocusing nonlinearity, respectively. Stationary solu- 
tions are found in the form En = Un{x) exp(icrfc„z) where 
kn are the propagation constants of the components. We 
assume strong saturation regime, s = 1, which is closer 
to realistic experimental conditions. 

The induced waveguiding regime, well studied in the 
context of vector solitons corresponds to the case 
when the intensities of the two interacting fields are sig- 
nificantly different. Then the strong field (e.g., ui) is de- 
scribed by a single (scalar) NLS equation, and the weaker 
field propagates in the effective linear waveguide induced 
by the stronger component via XPM. Here we assume 
that the effective waveguide (i.e. the grating) is created 
by a periodic nonlinear field ui{x) with the propagation 
constant fci , described by the stationary wave-train solu- 
tions of a scalar equation ^ (see also Ref. |^). Integrat- 
ing the stationary form of Eq. |^ once, we introduce the 
effective potential P{ui) = (cr/2)[(l — fci)?/^ — ln(l -t-u^)], 
so that the general stationary solution ui {x) with the am- 
plitude A can be found by solving the equation P{A) = 
^{dui/dx)^ + P{ui). Figure^a) shows the form of the 
potential P{ui) in both focusing (cr = +1, top) and de- 
focusing (cr ~ —1, bottom) cases. The minima of P{ui) 
correspond to a plane wave with the constant amplitude 
A^^ = ki/{l — fci), whereas the bright soliton solutions 
correspond to the separatrix at A = As with P{As) = 0. 
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FIG. 1: (a) Effective potential P(u\) for the focusing (top) 
and defocusing (bottom) cases (k\ = 0.5). Period (b) and 
examples (c-e) of the refractive index modulation A'^(/) for 
the three branches of periodic solutions of the NLS equation: 
(c) = 0.8, (d) = 1.41, and (e) A2 = 2.01. 
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FIG. 2: The bandgap structure of the linear spectrum 
—k2{A), induced by the nonlinear periodic grating ui{x) for 
ki = 0.5. The bands are shaded. Marked dots indicate the 
propagation constants for gap solitons shown in Fig. |3 



In the limit s ^ 0, exact analytical expressions for the 
nonlinear periodic waves ui(x) can be written down for 
a = ±1, in terms of the elliptic Jacobi functions, cn(x, /i), 
dn(a;,/i), and sn(a:,/j,), with the modulus < ^(fci) < 1. 
They represent self-consistent solutions of the cubic NLS 
equations which coincide with well-studied Hill's equa- 
tion with associated Lame potentials 0. It can be 
shown that the general structure of the periodic solu- 
tions is preserved for s 0. For a ~ +1, there exist 
two branches of the periodic (cnoidal) solutions shown 
through their induced refractive index modulation in 
Fig.IHd), for A = Ai, A^^ < Ai < As, and in Fig. [He), 
for A = A2 > As- The cn-type solutions of the branch A2 
have nodes, whereas the dn-type solutions of the branch 
Ai are nodeless. In the defocusing case {a = —1), there 
exists only one branch of the srt-type periodic solutions 
for A ^ Aq < Acv,, see Fig. ^c). In a strongly nonhn- 
ear limit the large-period Ai^2 and ^0 solutions describe 
periodic trains of bright and dark solitons, respectively. 



Having identified the stationary, nonlinear periodic 
solutions for the scalar field ui{x), we find that, in 
the induced waveguiding regime, the weak wave U2 is 
scattered by an effectively fixed linear grating charac- 
terized by the potential N{I), where / = uf(a;). The 
guiding properties of such a linear grating are deter- 
mined by the bandgap structure of the spectrum of 
the Hill's equation: (fu2/dx'^ = ~aN{T)u2 + ^2^2, 
where the eigenvalue k2{A) depends on the grating am- 
plitude. The eigenfunctions satisfy the Bloch condition 
U2ix) = BKpiiK L)u2{x + L), where L is the period and 
K is the momentum of the lattice. The spectrum con- 
sists of M bands and a continuum band, with the total 
TO = 2M -I- 1 band edges. The eigenfunction at the m-th 
band edge corresponds to a strictly periodic Bloch wave 
= fe"(a;) = ±b"'{x + L), for which KL = 0,7r. Fig- 
ure [3 shows an example of the bandgap spectrum —k2 
generated by the scalar cnoidal wave ui{x), for ki — 0.5. 
The Bloch waves at the band edges to = 1, 2, 3, 4, 5, . . . 
have the propagation constants —ak^ff ^ "O-^™^^ and 
the period L, 2L, 2L, L, L, . . .. The Bloch wave b"^{x) at 
the band edge fc™ = ki coincides with the scalar cnoidal 
wave Ml (x) . In the case of a saturable nonlinearity some 
predictions of the number and position of the bands and 
gaps can be made using theory of Lame-type equations 
[31 . For example, in the case of a defocusing nonlinear- 
ity, the grating potential can be well approximated by 
A''(a;) « a{a + l)/isn^(a;, /i), where a = 1. Since the spec- 
trum of the Lame's equation has M — a bound bands, 
it is expected that the Ag-tyjie grating generates single 
bound band followed by a semi-infinite band (as seen in 
Fig.[2|for A<^cw). 

In the limit A ^ As the period of the cnoidal-wave 
solution diverges and the spectrum bands disappear, see 
Fig. [2 From the other hand, when the grating ampli- 
tude A approaches the plane-wave amplitude Acv, in the 
focusing case, the periodic modulation of the refractive 
index vanishes and the gaps disappear, see Fig. [21 

To be useful for creation of robust dynamical photonic 
structures, nonlinear periodic waves should be stable. 
Previous studies of stability of periodic solutions 0, S 
suggest that the solutions of the Ai-type are strongly 
unstable due to modulational instability (MI), whereas 
MI is suppressed for the A2-type solutions in a saturable 
medium, and also for Ao-type solutions in the defocus- 
ing case. Our numerical studies have confirmed that the 
^Q-type grating in the defocusing case is both linearly 
and dynamically stable, and also demonstrated that the 
^2-type solutions are only weakly (oscillatory) unstable. 
In contrast, the Ai-type lattice is quickly destroyed by 
strong symmetry-breaking instabilities; therefore we ex- 
cluded it from our further consideration. 

The localization of the probe field U2 in the gaps of the 
linear spectrum of the periodic structure induced by the 
field ui{x) can occur in the nonZmear regime of the probe 
propagation through the grating. In this regime, the sig- 
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FIG. 3: Examples of the stationary two-component solutions 
in the different bandgaps for fixed parameters of the periodical 
component ki = 0.5 for the defocusing medium (the upper 
row) and the focusing medium (two lower rows). Dashed - 
the periodic component ui; solid - the localized component 
U2 ■ Gap solitons correspond to the marked points in Fig. 2] 



nificant intensity of the probe beam does not permit to 
neglect its nonlinear self-action. When the back-action 
of the probe on the grating through XPM is ignored (e.g. 
in the case of a weak material nonlinearity for the grat- 
ing wave) , the physics of the localization is similar to the 
standard case of nonlinear waves in fixed periodic po- 
tentials, well studied in the context of both optical and 
matter waves 0, H, 0, 0. However, in our problem 
the grating and scattered wave are strongly nonlinearly 
coupled and, therefore, as in the case of two-component 
vector solitons 0] , we should expect the existence of self- 
consistent hybrid structures formed by a periodic wave 
and a localized gap mode. 

Indeed, by solving vector Eq. numerically, with a 
value of fci fixed to that of the scalar grating ki = kf, 
we have found different families of solutions of Eq. 
consisting of the oscillatory (ui) and localized (U2) mu- 
tually trapped components. The propagation constant 
of the localized component always lies within the gaps of 
the linear spectrum. Therefore this component can be 
described as a gap soliton with even or odd symmetry 
centered at a maximum or minimum of the grating 
potential, respectively. Figure 13 shows some of examples 
of such a gap solitons for both defocusing and focusing 
cases. First, we note that the powers of discrete solitons 
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FIG. 4: Power of the localized component, Q = J U2{x)dx, 
vs. k2 for ki — 0.5 and (top) defocusing (^0 = 0.8) and 
(bottom) focusing {A2 — 2.01) nonlinearity. Shaded regions 
corresponds to the spectral bands, and the dots mark the 
solutions shown in Fig. |21 



with different symmetries coincide, i.e. these solitons be- 
long to the same family. This indicates the absence of 
the Peierls-Nabarro potential barrier and good mobility 
of the localized states. Second, due to the nonlinear XPM 
interaction, the induced grating is strongly modified by 
the localized component, but recovers periodicity in the 
far field. Different cases of the gap solitons are summa- 
rized in Fig. 01 where we show the families of localized 
modes for both focusing and defocusing cases. 

In agreement with the theory of gap solitons in non- 
linear periodic structures the families of localized 
states originate at the edges of the bands with the num- 
bers TO = 2 (for cr < 0) or m = 1, 3, 5... (for cr > 0), where 
the effective dispersion, {d'^k2/dK^)\k^^ , is correspond- 
ingly negative or positive. At the respective band edge, 
the low-power gap soliton is weakly localized, and can be 
described as a slowly varying envelope of the correspond- 
ing Bloch wave b"^{x) [iflj. In the defocusing case, only 
gap solitons [(a) and (b)] can exist in the induced grating, 
whereas in the focusing case, both gap [(c) and (d)] soli- 
tons and self-trapped [(e) and (f)] solitons in the semi- 
infinite gap are possible. Near the opposite gap edges, 
where the gap modes have high powers, the periodic wave 
of the grating acquires significant defects induced by the 
localized state, however both components still exist as a 
vector stationary state. 

To understand the nature of this composite state, we 
consider the correction to the linear grating spectrum 
due to the low-amplitude gap mode in the U2 compo- 
nent, which is bifurcating off the lower edge of the band 
II (to — 3) in the focusing case, or upper edge of the 
band I (to = 2) in the defocusing case. Near the bi- 
furcation threshold, the nonlinear XPM coupling leads 
to the effective shift of the propagation constant of the 
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FIG. 5: Stable propagation of the gap soliton in the defocus- 
ing grating, (a) Propagation dynamics of the initial (odd) 
state, corresponding to the family (a-b) in Fig. 2] perturbed 
by a random amplitude noise, in the absence of the grating; 

(b) the gap mode (solid) and grating (dashed) initial profiles; 

(c) propagation dynamics in the presence of the grating; (d) 
the final state at z = 400. 



firmed dynamical stability of band-gap solitons by nu- 
merical integration of the vector dynamical model 1^. 
Figure shows an example of the stationary propaga- 
tion of an odd gap soliton [family (a-b) in Fig. 0] in an 
induced nonlinear grating for the defocusing case. Both 
components are initially perturbed by a random noise at 
20% of their peak amplitude [Fig. O (b)]. If the grat- 
ing is removed, the localized gap mode can no longer be 
supported by the defocusing nonlinearity and strongly 
diffracts [Fig. (a)]. Being coupled to the lattice, the 
gap mode generates a defect in the grating and coexists 
with it as a dynamically stable composite state, which is 
clearly robust to perturbations [Figs. El (c,d)]. 

In conclusion, we have introduced novel composite 
band-gap solitons where one of the components creates a 
periodic nonlinear lattice which localizes the other com- 
ponent in the form of a gap soliton. Nonlinear local- 
ization of this kind should be generic to models of non- 
linear ly interacting multi-component fields, where one of 
the components can exist in a dynamically stable self- 
modulated periodic state. Here, we considered a spe- 
cific example of a spatial nonlinear photonic structure 
induced by optical beams in a photorefractive crystal. 
Another example is the dynamical Bragg gratings for op- 
tical pulses obtained through the cross-phase modulation 
in highly birefringent fibers 
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periodic grating component: ki = kf + Aki, where fcf 
is the fixed propagation constant for the scalar grating 
Ui{x), and Afci ~ cr / u^u^l + u\)~'^dx. As a result, 
the band edges fc™ corresponding to the periodic Bloch 
modes ¥^{x) = u\{x) shift into the gap. This means 
that the nonlinear mode ui, coupled to the field U2 and 
corresponding to the fixed propagation constant fci — kf 
(i.e. fixed power), now lies within the band, and has an 
oscillatory, but not strictly periodic nature. Therefore, 
the localized low-power gap mode in the U2 component is 
nonlinearly coupled to an in-band dark-like defect mode 
in the ui component, which is localized on the Bloch- 
wave background. Together, the two components form a 
novel band-gap composite soliton. In the deeply nonlin- 
ear regime, when the high-intensity gap mode induces a 
large defect in the grating wave, both components exist 
as a vector soliton in a mutually formed periodic waveg- 
uide N{I), which now depends on both fields intensities, 
I = + u\. As ^2 — > ki, the gap mode's amplitude 
approaches that of the grating, and both components of 
the band-gap soliton become oscillatory. 

The crucial issue of stability of the gap solitons in the 
nonlinear induced gratings is therefore linked to the sta- 
bility of the composite band-gap states. We have con- 
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